Bose-Einstein condensate in a harmonic trap decorated with Dirac 6 functions 
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We study Bose-Einstein condensation in a harmonic trap with a dimple potential. We specifically 
consider the case of a tight and deep dimple potential which is modelled by a Dirac 8 function. This 
allows for simpler, explicit numerical and analytical investigations of noninteracting gases. Thus, 
the Schrodinger equation is used instead of the Gross-Pitaevski equation. Calculating the atomic 
density, chemical potential, critical temperature and condensate fraction, the role of the relative 
depth of the dimple potential with respect to the harmonic trap in large condensate formation at 
enhanced temperatures is clearly revealed. 

PACS numbers: 03.75.Hh, 03.65.Ge 



I. INTRODUCTION 

Bose-Einstein condensation was the last major discov- 
ery of Einstein [J, By applying the Bose-Einstein 
statistics to an ideal gas, Einstein showed that "from a 
certain temperature on, the molecules condense without 
attractive forces .. . " |3| and discovered the Bose-Einstein 
condensation in 1925. This theoretical prediction moti- 
vated the experimental studies for realizations of Bose- 
Einstein condensates of gases. Seventy years after the 
prediction of Einstein, Bose-Einstein condensates (BEC) 
of dilute gases have been observed at very low tempera- 
tures by using ingenious experimental designs [4|, [5j, l6( . 

Experimentally available condensates are systems with 
finite number of atoms N, confined in spatially inhomoge- 
neous trapping potentials. Their understanding requires 
theories that go beyond the usual treatments based upon 
London's continuous spectrum approximation [7j 
thermodynamic limit N — > oo. Such studies 
veal that BEC can occur in harmonically trapped lower- 
dimensional systems for finite N despite the enhanced 
importance of phase fluctuations [l(J . Quasicondensates 
with large phase fluctuations ma y st ill occur [ll|. This 
is in contrast to standard results [12| , in agreement with 
the Mermin-Wagner-Hohenberg theorem [l3|> [l4| in the 
thermodynamic limit. Bose-Einstein condensation in one 
dimension (ID) with a harmonic trap is attractive due to 
enhanced critical temperature and condensate fraction 
[1, Q . Recently, one- and two- dimensional BECs have 
been created in experiments [Hj]. One-dimensional con- 
densates were also generated on a microchip [lj| [TtJ and 
in lithium mixtures [18j ]. 

Modification of the shape of the trapping potential 
can be used to increase the phase space density 
"Dimple" -type potentials are the most favorable poten- 
tials for this purpose [2(| Hi], H3j- The phase-space 
density can be enhanced by an arbitrary factor by us- 
ing a small dimple potential at the equilibrium point of 
the harmonic trapping potential [20]. A recent demon- 
stration of caesium BEC exploits a tight dimple poten- 



tial [2l|. Quite recently, such potentials were proposed 
for efficient loading and fast evaporative cooling to pro- 
duce large BECs [23|. Tight dimple potentials for one- 
dimensional (or strictly speaking quasi-one-dimensional) 
BECs offer attractive applications, such as controlling 
interaction between dark solitons and sound [24j . intro- 
ducing defects such as atomic quantum dots in optical 
lattices [25} . or quantum tweezers for atoms 26]. Such 
systems can also be used for spatially selective loading of 
optical lattices 27]. In combination with the condensates 
on atom chips, tight and deep dimple potentials can lead 
to rich novel dynamics for potential applications in atom 
lasers, atom interferometers and in quantum computa- 
tions (see Ref. 



281 ] and references therein) 



In this paper, we investigate Bose-Einstein condensa- 
tion in a ID harmonic trap decorated by a tight and 
deep dimple potential. The quantum kinetics of a sim- 
ilar system but for a 3D single-mode trap, modelled by 
a deep but narrow spherical square potential well, was 
studied in Ref. [29j . where a master equation, to de- 
scribe condensate growth, was developed. We model 
the dimple potential using a Dirac 5 function. The 
delta function can be defined via a Gaussian function 
[30T ] g(x,a) = (I/^/tto) cxp (— x 2 /a 2 ) of infinitely narrow 
width a so that g(x, a) — > 6(x) for a — > 0. This allows 
for analytical calculations in some limiting cases as well 
as a simpler numerical treatment for arbitrary parame- 
ters. We calculate the transition temperature as well as 
the chemical potential and condensate fraction for vari- 
ous number of atoms and for various relative depths of 
the dimple potential. For describing a system with in- 
teracting particles, the Gross-Pitaevski equation is usu- 
ally utilized. We note that we neglect the interactions 
between atoms in our model, and thus the Schrodinger 
equation for the harmonic trap with the dimple potential 
is solved. 

The paper is organized as follows. In Sect. II, we 
present the analytical solutions of the Schrodinger equa- 
tion for a Dirac (5-decorated harmonic potential and the 
corresponding eigenvalue equation. In Sect. Ill, deter- 
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mining the eigenvalues numerically, we show the effect of 
the dimple potential on the condensate fraction and the 
transition temperature. In Sect. IV, we present a semi- 
classical method to calculate condensate fraction when a 
dimple potential is added to harmonic trap adiabatically. 
Analytical results in the limit of strong dimple potential 
are presented in the Sect. V. Finally, we present our 
conclusions in Sect. VI. 



II. HARMONIC POTENTIAL DECORATED 
WITH DIRAC DELTA FUNCTIONS 

We begin our discussion with the one dimensional har- 
monic potential decorated with the Dirac 8 functions [3lj ] - 
[33j |. This potential is given as 
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where uj is the frequency of the harmonic trap, P is a 
finite integer, and oVs are the strengths (depths) of the 
dimple potentials located at Xi's with xi < X2 < ■■■ < 
xp with Xi 6 (— oo,oo). The factor ti 2 /2m is used for 
calculational convenience. A negative Oi value represents 
repulsive interaction while positive <7j value represents 
attractive interaction. We can write time-independent 
Schrodinger equation as 
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+ V(x)^{x) = EW(x) 
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By inserting E = (£ + ^)huj, with £ a real number, and 
introducing dimcnsionlcss quantities z = x/xq, and Zi = 
Xi/xo with xq = y/h/2rnu>, the natural length scale of 
the harmonic trap, we can re-express Eq. ([2]) as 
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where Aj = x^Ui. For z ^ Zi, the Eq. ((2|) has two lin- 
early independent solutions. For £ ^ 0, 1, 2, ... these lin- 
early independent solutions are parabolic cylinder func- 
tions D^{z) and D^(—z), denned as: 
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where $ (-£/2, 1/2; z 2 /2) and $ ((1 - £)/2, 3/2; z 2 /2) 
are the confluent hypergeometric functions $ (a, b; z 2 /2) 



defined in Ref. [3J for a = £/2, b = 1/2 and a = 
(1 — £)/2, & = 3/2, respectively. Thus, the general so- 
lution of Eq. Q between the locations of two Dirac 
8 functions is 



= a t D ( (z) + b l D ( {-z), 
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with di and bi can be written in terms of b\ via transfer 
matrix method and b\ is determined by normalization 
(ai = fep+i = due to square integrability) [H, 
D^(—z) is regular as z — > — oo, but |-Df(— z)\ — ► +oo 
as z — > oo, and D^(z) is regular as z ^ oo, but 
|Dj(z)| — * +oo as z — > — oo [3J|. By using the trans- 
fer matrix method explained in Refs. (3ll. l35j|. one finds 
the eigenvalue equation for the potential given in Eq. |T]) . 
For the special case of the potential given in Eq. ([I]) , har- 
monic potential decorated with one Dirac 8 function at 
its center P = 1, x\ = 0, 
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the eigenvalue equation of the even states is written as 



T((l-0/2) _ r(3/4- E/2hw) _ ay/h/muj 
r(-£/2) ~ T(l/4-E/2ha;) ~ 4 ' 



(8) 



In this case, the odd solutions of the harmonic potential 
are unaffected because the Dirac 8 function is located 
at the center of the harmonic potential. On the other 
hand, the energy eigenvalues of even states change as 
a function of a. The ground state energy eigenvalue de- 
creases unlimitedly as a increases (attractive case). How- 
ever, the energies of the excited even states are limited 
by the energies E^n+i of odd states and as er — > oo, 
E2n+2 — * E 2n +i = {^n + 1 + l/2)hjj where n = 0, 1, .... 
Thus, as a — > oo the energy eigenvalues of these states 
go to the energy eigenvalue of the one lower eigenstate, 
and the odd energy eigenstates asymptotically become 
doubly degenerate. 



III. BEC IN A ONE-DIMENSIONAL 
HARMONIC POTENTIAL WITH A DIRAC 
8 FUNCTION 

We begin our discussion about BEC in a one- 
dimensional harmonic potential decorated with a Dirac 
8 function by investigating the change of the critical tem- 
perature as a function of a. One estimates a er value using 
the parameters of Ref. [22| . In this paper, the minimum 
value of the dimple potential is given as U c — ks 4 /iK 
and average potential width is given as r = 1 — 100 /im. 
We equate the strength (—H 2 a/2m) to the product U c r 
to get an estimate of a. We find that a varies approxi- 
mately between 10 s 1/m and 10 10 1/m as r changes from 
1 /imto 100 /im. We define a dimensionless parameter in 
terms of a as: 
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FIG. 1: The critical temperature T c and T c sc vs A for AT = 
10 4 . A is a dimensionless variable defined in Eq. (JSJ. The 
crosses and dots show T c and T c sc values, respectively. Here 
we use m = 23 amu ( 23 Na) and u = 2vr x 21 Hz [l3|. The 
logarithmic scale is used for A axis. 



If 10 8 l/m < a < 10 10 l/m then 460 < A < 46000 
for the experimental parameters m= 23 amu ( 23 Na), 
u = 2 vr x 21 Hz [33] and 230 < A < 23000 for the exper- 
imental parameters m= 133 amu ( 133 Cs), ui = 2tt x 14 
Hz [21] . In this work, we show that, even for small A 
values, the condensate fraction and critical temperature 
change considerably. 

The critical temperature (T c ) is obtained by taking 
the chemical potential equal to the ground-state energy 
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at T = T c , where /3 C = l/(k B T c ). For finite N value, we 
define T c ° as the solution of Eq. (ID]) for A = (only the 
harmonic trap). 

In Eq. (fTOj) . Ej's are the eigenvalues for the potential 
given in Eq. ([7]). The energies of odd states are un- 
changed and equal to (2n + 1 + \j2)Hui. The energies 
of even states are found by solving Eq. ([5]) numerically. 
Then, these values are substituted into the Eq. (fTU)) . and 
finally this equation is solved numerically to find T c . We 
obtain T c for different A and show our results in Fig. [T] 
In this figure, logarithmic scale is used for A axis. As A 
increases, the critical temperature increases very rapidly 
when A > 1. The solid line in Fig. Q] shows the change of 
the critical temperature with A. Here we take N = 10 4 
and use the experimental parameters m= 23 amu ( 23 Na) 
and uj = 2?r x 21 Hz [37j . 

We have calculated the change in the critical tem- 
perature approximating dimple-type potentials by Dirac 
S function. The change in the critical temperature can 



FIG. 2: (T c sc -T c /r c ) x 100 vs A, the difference of the critical 
temperatures calculated by the semiclassical and our method 
with respect to A. The logarithmic scale is used for A axis. 



also be calculated by assuming that the optical trap is 
added to harmonic trap adiabatically [2fj, [38|. In this 
semi-classical approach, one uses the density of states 
and assumes the energy spectrum is continuous utilizing 
the fact that fiw/^T) <C 1. For a one dimensional har- 
monic trap, the number of particles in the thermal gas 
can be found as 



N = N + 
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where iVo denotes the average number of particles in the 
ground state. The 1 / (Huj) factor before the integral in Eq. 
(|11[) is the density of states for one dimensional harmonic 
potential. The critical temperature T~ in the harmonic 
trap can be found by taking = Huj/ 2 and writing the 
following equation for the total number of particles, N: 
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where 0® = 1/ (ksT®) . The addition of an attractive dim- 
ple potential to the harmonic trap decreases the ground- 
state energy, resulting a decrease in the chemical poten- 
tial. Since the chemical potential is assumed to be equal 
to the ground sate energy at critical temperature, replac- 
ing fi in Eq. (fTT|) by E g , the ground state energy for the 
harmonic trap together with the dimple potential, it is 
possible to calculate the critical temperature of the new 
system; 
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where j3^ c = l/(fcsT ( f c ') denotes the critical tempera- 
ture values obtained by semiclassical method. Here we 
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FIG. 3: The chemical potential fj, vs temperature T /T® for 
TV = 10 4 and solid line for A = and dashed line for A = 46. 
The other parameters are the same as Fig. [1] 
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FIG. 4: N /N vs T/T c ° for TV = 10 6 and solid line for A = 0, 
dashed line for A = 4.6, dotted line for A = 46. The other 
parameters are the same as in Fig. [1] 



assume that the energies of the excited state are same 
with the harmonic trap by following Ref. ^20]. We have 
calculated the critical temperature for various dimple po- 
tentials differing according to potential depths. We take 
the Eg values as the ground state of the potential given 
in Eq. (|7|) for different a [ or A = a^Jhjirnio) ] values 
in order to compare the results of this method with our 
results. 

The critical temperature values calculated using this 
semiclassical approach (T^ ) is shown in Fig. Q]as points. 
As seen from the figure, the critical temperature val- 
ues estimated by our method and by the semiclassical 
method are almost equal. We show in Fig. {2} percent- 
age difference of the critical temperatures that are calcu- 
lated by using these methods. The difference of these two 
methods for very small A values is due to the accurate 
calculation of the summation with accurate energy val- 
ues for the excited states for our method; however, the 
summation is approximated with an integral by using 
the harmonic trap density of states for the semiclassical 
method. 

For a gas of N identical bosons, the chemical potential 
fi is obtained by solving 

N = Y~-51 S = ^0 + V-57 ; , (14) 

i—0 i—1 

at constant temperature and for given N, where e% is the 
energy of state i. We present the change of /i as a function 
of T/T c ° in FigEJfor TV = 10 4 ; A = and A = 46. By 
inserting fi values into the equation 



we find the average number of particle in the ground 
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FIG. 5: Condensate N /N vs temperature T/T° for A = 46. 
The solid line for N = 10 s , dashed line for TV = 10 6 , dotted 
line for TV — 10 4 . The other parameters are the same as in 
Fig. ID 



state. TVo/TV versus T/T° for TV = 10 6 and A = 0,4.6,46 
are shown in Fig. |H In this figure, the result for A = 
is the same as the result obtained by Ketterle and van 
Druten 0]. As mentioned in Ref. Q, the phase transi- 
tions due to discontinuity in an observable macro param- 
eter occurs only in thermodynamic limit, where TV — * oo. 
However, we make our calculations for a realistic system 
with a finite number of particles. Thus, TV /TV is a finite 
nonzero quantity for T < T c without having any discon- 
tinuity at T = T c . 

It is useful to know the behavior of the condensate frac- 
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FIG. 6: Condensate fraction No/N vs the strength of the 
Dirac S potential A when N = 10 4 . The logarithmic scale is 
used for the A axis. The other parameters are the same as in 

Fig.LU 



tion as a function of temperature for a fixed value of A. 
We present the condensate fraction for N — 10 4 , 10 6 , 10 8 
when A = 46 in Fig. All condensate fractions for dif- 
ferent N values are drawn by using their corresponding 
T c ° values. These T c ° values are 13 fiK for N = 10 4 , 85 fiK 
for N = 10 6 and 6200 for N = 10 8 . 

We also find the condensate fraction as a function of 
A at a constant T = T c °. These results for N = 10 4 are 
shown in Fig. [Sj In the following section we show that 
the condensate fraction changes exponentially for large 
A. Thus, large A values (A > 1) induce sharp increase in 
the condensate fraction. 

Finally, we compare the density profiles of condensates 
for a harmonic trap and a harmonic trap decorated with 
a S function (A = 4.6) in Fig. [7] Since the ground-state 
wave functions can be calculated analytically for both 
cases, we find the density profiles by taking the absolute 
square of the ground-state wave functions. 



IV. ADIABATIC CONDENSATION 

One can also find the condensate fraction assuming 
a dimple potential is added to harmonic trap adiabati- 
cally [20, [38|] . We first calculate the grand potential of a 
one-dimensional Bose-gas in a harmonic trap using semi- 
classical approximation: 
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where O denotes the contribution of the ground state 
to the grand potential. The l/(fko) factor before the 
integral in Eq. (fH?)) is the density of states for one- 



FIG. 7: Comparison of density profiles of a BEC in a har- 
monic trap with a BEC in a harmonic trap decorated with a 
S function (A = 4.6). The solid curve is the density profile 
of the BEC in decorated potential. The dashed curve is the 
density profile of the ID harmonic trap (A = 0) . The param- 
eter z is dimensionless length defined after Eq. @. The other 
parameters are the same as in Fig. [1] 



dimensional harmonic potential. We calculate the en- 
tropy using S = —dQ/dT and get 
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where 52 ( x ) is a Bose function and is defined as 
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For temperatures slightly above Tj, the entropy is 
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assuming /j — fko/ 2. We assume that adding a dim- 
ple potential changes the ground state and the chemical 
potential such that // = E g < hu/2 and the entropy ex- 
pression becomes 



j5 

k B 



2kjiT g* (c 0(E ^ 
fko 



£-1"" 



1 -e 



/3(E g -3huj/2) 



(20) 



Since the entropy remains constant during an adiabatic 
process, we can use Eqs. (|19p and (|20p to calculate the 
final temperature of the system. Then, we calculate 
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FIG. 8: The comparison of the condensate fraction values 
evaluated with the semiclassical method and our method for 
N — 10 4 . The crosses and dots show the results of the semi- 
classical method and our method, respectively. The logarith- 
mic scale is used for the x axis, where AE = hu/2 — E g . 



the condensate fraction for this final temperature. We 
present the change of condensate fraction with respect 
to AE/(k B T°) in Fig. E] where AE = hj/2 - E g . The 
solid and dashed lines in Fig. [5] show the results of the 
semiclassical method and our method, respectively. AE 
increases with A, representing the potential depth for a 
constant width. We have calculated the condensate frac- 
tion for the same ground-state energy values of the dim- 
ple potential in both the semiclassical and our method 
to make a comparison. The condensate fraction results 
in Fig. [8] differ for these two methods for very small A 
values. This is due to two effects: (a) Calculations of the 
chemical potential is done accurately by using the con- 
stant N for the finite number of particles for our method, 
however the chemical potential is just taken the ground- 
state energy for the semiclassical method 0, H^]. (b) 
An accurate summation calculation is performed for our 
method; however, the summation is approximated by an 
integral for the semiclassical method. 

The plots in Figs. Q] and [8] show that the agreement of 
the results of two methods is better for large A (or large 
AE) values. This is due to the fact that the average 
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occupations for excited states is so low (n 
that the semiclassical treatment for this gas is valid. 



V. APPROXIMATE SOLUTIONS OF THE 

CRITICAL TEMPERATURE AND 
CONDENSATE FRACTION FOR LARGE a 

If we apply a deep dimple potential (p — > oo) to the 
atomic condensate in a harmonic trap, the problem can 
be solved analytically by approximating the summation 



for N with an integral. As T — > T c , /i 



Eg, where 
(21) 



and it is the bound state energy of a Dirac 6 poten- 
tial HI]. By using density of states p(e) and utilizing 
E2n+2 — * E^n+i (odd-state energies) as a — > oo, the 
summation in Eq. (|14p is converted to the following in- 
tegral 
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calculating the integral, we get 
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Defining A — \Es\/Huj, we have A ^> 
exp ((— /3 c hujA)) 1 for very large a. Then, 
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for this case and we get for the critical temperature 
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for large A/N value. For one-dimensional experimental 
systems, N w 10 3 — 10 6 atoms. For a specific case A/N « 



10 4 
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which shows that the critical temperature increases lin- 
early with increasing bound-state energy for the dimple 
potential. 

By using Eq. (|24[) . the condensate fraction can be 
written as 



iV ° = i J_ e -l3f3 c k B (T c -T)huA 
N T r . 



(27) 



for T < T c . This result indicates an exponential increase 
of Nq as a function of T for T < T c . Thus, the number 
of atoms in Bose-Einstein condensate will rise drastically 
when a very strong, very short-range (point) interaction 
is added to a harmonic confining potential. However, we 
ignore the interactions between the atoms and neglect 
nonlinear terms in the Gross-Pitaevski equation which 
will modify these results. 



VI. CONCLUSION 

We have investigated the effect of the tight dimple 
potential on the harmonically confined one-dimensional 
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BEC. We model the dimple potential with the Dirac 
8 function. This allows for analytical expressions for 
the eigenfunctions of the system and a simple eigenvalue 
equation greatly simplifying numerical treatment. Pure 
analytical results are obtained in the limit of infinitely 
deep dimple potential case. 

We have calculated the critical temperature, chemi- 
cal potential and condensate fraction and demonstrated 
the effect of the dimple potential. We have found that 
the critical temperature can be enhanced by an order of 
magnitude for experimentally accessible dimple potential 
parameters. In general, T c increases with the relative 
strength of the dimple potential with respect to the har- 
monic trap. In our model system, the increase in the 
strength of the Dirac 8 function can be interpreted as 
increasing the depth of a dimple potential. 

In addition, the change of the condensate fraction with 
respect to the strength of the Dirac 8 function has been 
analyzed at a constant temperature (T = T°) and with 
respect to temperature at a constant strength. It has 
been shown that the condensate fraction can be increased 
considerably and large condensates can be achieved at 
higher temperatures due to the strong localization effect 
of the dimple potential. Analytical expressions are given 
to clarify the relation of the condensate fraction and criti- 
cal temperature to the strength of deep dimple potential. 

We also show that similar results for the critical tem- 
perature and condensate fraction can be obtained by us- 



ing a semiclassical approximation. However, in this sim- 
pler approach, it is not possible to determine the density 
profile and when the thermodynamic limit is not satis- 
fied, the estimations of semiclassical approach fail. 

Finally, we have determined and compared the density 
profiles of the harmonic trap and the decorated trap with 
the Dirac 8 function at the equilibrium point using an- 
alytical solutions of the model system. Comparing the 
graphics of density profiles, we see that a dimple poten- 
tial maintains a considerably higher density at the center 
of the harmonic trap. 

The results presented are obtained for the case of non- 
interacting condensates for simplicity. The treatment 
should be extended to the case of interacting conden- 
sates in order to make the results more relevant to experi- 
mental investigations. This case deserves further detailed 
and separate calculations. We believe the method pre- 
sented of (5-function modelling of tight dimple potentials 
can help significantly such theoretical examinations. 
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